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' A generalizing formulation of dynamical real-space renormalization that suits for arbitrary spin 

' systems is suggested. The new version replaces the single-spin flipping Glauber dynamics with the 

single-spin transition dynamics. As an application, in this paper we mainly investigate the critical 
slowing down of the Gaussian spin model on three fractal lattices, including nonbranching, branching 
and multibranching Koch Curve. The dynamical critical exponent z is calculated for these lattices 
using an exact decimation renormalization transformation in the assumption of the magnetic-like 
perturbation, and a universal result z = 1/u \s found. 
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I. INTRODUCTION 



' The dynamics of spin systems approaching their second-order phase transition points have been an important subject 
of many studies in the last few decades. One of the interesting phenomena is the critical slowing down characterized 
i' by a divergent relaxation time r. A reasonable explanation is seemingly that the long raiige fluctuation leads to long 
'"^ . time evolution of the order parameter. According to the dynamical scaling hypothesis the divergent relaxation 
^ ' time T and the divergent correlation length C, can be related by r ~ , where z is the dynamical critical exponent and 
O ' is believed to depend only on large universal features of the model Hamiltonian and the assumed dynamic process . 

Obtaining the exact solution based on a master equation, except for a few cases [^-^, is not an easy job. One have 
to evaluate it by means of approximate methods, such as the Monte Carlo simulation, the high-temperature series 
expansion, etc. However, the success of renormalization-group (RG) methods in obtaining the critical exponents 
^ ■ and universality classes of static problems led to several attempts to use RG ideas in critical dynamics. One of the 
typical examples is a generalization of the e expansion technique in which the dynamics is described by a Langevin 
equation. This method enables the calculation of the time-dependent correlation functions, but it is only useful near 
the upper critical dimension. Another Q is a generalization of the real-space RG techniques, the starting point is 
a master equation instead of Hamiltonian. This method is preferable in discrete spin systems, and could be used to 
fSl , directly calculate the dynamical critical exponent. In addition, it is simple and transparent, and very accurate for 
. certain systems so that it has been used quite extensively in the last years ]9|-p^. Other examples see Refs. |16 1^. 



The dynamical real-space renormalization group (DRSRG) technique, proposed by Y. Achiam and J. M. Kosterlitz 



and perfected by D.Kandel |18|, is our focus in this paper. First, we establish a formulation of DRSRG applying 
Q ' to arbitrary spin systems. Then, we investigate the critical slowing down of the continuous spin model on different 
r fractal lattices. In the generalizing formulation of DRSRG, we replace the single-spin flipping Glauber dynamics |^ 
^ ■ with the single-spin transition dynamics , and use the same notation of Rcf. [|l8| to express the critical dynamical 
Q ' exponent z. 

O ■ During the last many years, scientific journals have published many papers concerning critical dynamics of discrete 
^ spin systems, but a systematic study of the critical dynamics in continuous spin systems is very lacked indeed. The 
purpose of our latest papers |^,^ and this work is a attempt filling this gap. This is just our main motivation. We 
realize the fact that, though the Gaussian model is certainly an idealization, it is interesting and simple enough 
to obtain some fundamental knowledge of dynamical process in cooperative systems. So this is a ideal dynamical 
model that interests us greatly. We also realize that, as an extension of Ising model, the Gaussian model shows 
much differences from Ising model in the properties of static phase transition, and yet its knowledge of the dynamical 
behavior is unclear. Within the framework of single-spin transition critical dynamics in our previous paper we 
have obtained the dynamical critical exponent of the Gaussian model, z = l/i/ = 2, at the critical point Kc ~ h/2d 
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based on rigorous calculation. This means that the dynamical exponent is highly universal on translational symmetric 
lattices. However, what is the dynamical exponent on dilational symmetric lattice systems? All of these motivate us 
to finish this work. 

This paper is organized as follows: Section ^ is a detailed description of the dynamical real-space renormalization 
group (DRSRG) technique in which the dynamics is described by a Markov process with the single-spin transition 
instead of the single-spin flipping. In Sec. the critical dynamics of Gaussian model on three fractal lattices 
is studied. We take the exact decimation transformation and calculate the dynamical critical exponent z in the 
assumption of the magnetic-like perturbation. Section IV is our summary and discussion. 



II. DESCRIPTION OF THE METHOD 



In the single-spin transition critical dynamics M , the master equation can be written as 



-P ({a} , t) = - 5] 5] (1 - W^. (a, ^ <7,) P ({a} , t) 



dt 

i iy 

where pi is the transition operation defined by 

Vi! (o'i,o-2,...,cri,...,crAr,i) = J [ax,02, ■ ■ ■ -.Oi, . . . .aj^ ,t) , 

and Wi {(7i —^ cfi) is the single-spin transition probability which satisfies the following restraint conditions: 
(a) Ergodicity: 



(2.1) 



(b) Positivity: 

(c) Normalization: 



Vcrj,<Tj : Wj{<Tj -^d-j)^ 0; 



(2.2(a)) 
(2.2(b)) 
(2.2(c)) 



(d) Detailed balance: 



in which 



Vcrj,CTj : 



Wjiaj Oj) Peq{cri, . . . , CTj, . . . , ctat) ' 



(2.2(d)) 



Peq {{t}) = \ exph/37^({a})], Z = Y, exp[-/3W({a})], 

{^} 

where P^q is the equilibrium Boltzmann distribution function, Z the partition function and 7i({cr}) the system 
Hamiltonian. A well-chosen form of the transition probability is 



1 



■ exp 



exp 








E<T. exp 





(2.3) 



In order to study the critical slowing down, we can limit ourselves to the relaxation of an infinitely small perturbation 
from equilibrium. Following Y. Achiam's idea , two selections can be considered, that is magnetic-like perturbation 



F(M,i)^ 



Peq {{^}) 



(2.4) 



or energy-like perturbation 
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Pea ({^}) , 



(2.5) 



where qi distinguishes between points which have different R (the order of ramification), (i,j) denotes a sum over 

nearest-neighbor pairs, and hq. and h ^. ar e th e re duced external fields. 

Based on these two considerations (2.4) or ( |2.5D , the master equation (2.1) takes the following forms 

(2.6) 



- 5] /i^ (i) a.cT.Pe, ({a}) = - E E ~ ^^^^ P-'i ({^}) 



respectively. We can express further the Eqs. (2_^) and (2/7) as the unitized matrix formulation 

it) A {a) {k, {a}) = -h it) n (fc, a) P,, {k, {a}) , 



(2.7) 



(2.8) 



where, h (t) is a row matrix, A (cr) and (fc, a) are column matri ces. 

The critical dynamical behavior of the system described by Eq. (2.8), can be studied using the dynamical real-space 
renormalization group (DRSRG) technique. The DRSRG is composed of two stages. The first stage is the rescaling 
of the space by 



X —^ x' — Lx. 



(2.9) 



which is performed using a RG transformation (such as decimation or site-block transformation), where L is the 
length-rescaling factor. For example, in the case of the decimation transformation, the spins are divided into two 
groups {a} and {/Lt} under the control of a decimation operator T {p., a), then a trace over the {a} is performed. The 
process of decimation for a spin fimction / ({c}) can be demonstrated as 



^ {/ ( W)} = E ^ / ({'^i) - / (M) • 

It is certain that we need to rescale the interaction parameter k = J /kpT, and the spin /i, i.e. 

k ^ k' — R (fc) k, /i /i' = ^ (k) ^, 



(2.10) 



(2.11) 



so as to keep on an invariant form of the prob ability dist ribution, P^^ (fc', {/i'}), where ^ (fc) is the spin-rescaling factor. 
With the decimation transformation ( ^.IC ) , Eq. ( |2.8| ) takes the form 



-h {t) -Ra (k) • A' k') P4 (fc', {m'}) = -h {t) -Ro (k) • (/i', k') Pi^ {k', {fi'}) , 

where A' (/i', k') and n' (fj,' , k') retain the original form. Taking the monomark 

h' {t,k') = h(t)-RA (k). 



(2.12) 



which can be regarded as a RG transformation of the dynamic parameter h (t), Eq. (2.12) can be rewritten as 

|h'(t,A:')-A'(^',fc')P; ik'Af^'}) 

= -h' (t, k') ■ [R^i (k) . Rj, (k)] n' if,', k') P; [k\ {^^'}) . (2.13) 

The second stage of the DRSRG is the rescaling of the time by 

t-^t' ^ L'H, (2.14) 



which should be resulted in that Eq. (2.12) is restored to the invariant form of the master equation (|2.8|) 
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— h' {t\ k') -A' k') (fc', {^'}) = -h' {f, k!) n' (^', k') {k', {m'}) . 



(2.15) 



We might encounter two different cases in carrying out the Eq. (2.15). It can be reaUzed via the following analyses. 

First, for some homogeneous lattices with same coordination number {qi = q, hq. — h), Ra (k) and Rjj (k) are only 
1 X 1 mat rices. When system approaches its critical point kc, A = Ra ( fcc) = const, u = Rfj (fee) = const, then from 
Eq. ( 2.13| ) we can see that the invariant form of the master equation ( |2.8|) can be restored by preforming the time 
rescaling 



t 



t' = L-H = 



t 



X/uj' 



and from here we can further obtain the dynamical critical exponent 

_ ln(A/w) 



InL 



(2.16) 



(2.17) 



Second, for some inhomogeneous lattices with different coordination number, Ra {k) and Ro (fc) are m x m square 
matrices in which the order m of the matrices depends on the number of the parameter hq- . In this case we have to 
look for the invariant form at the limit of the order n — > oo of the RG transformation. Because our starting point 
is very close to the fixed point kc of the static RG transformation, the eigenvalues of t he tr ansformation matrices 
Ra (fc ~~* kc) and Rn (fc kc) control the scaling properties as i ^ oo ]To| , p^ . Hence, (2.17) again determines the 
dynamic exponent, merely X/uj should be replaced by Amax/wmin lH], i-e. 



In (Aniax/ "-^min) 
h^l^ 



(2.18) 



where A^ax is the largest eigenvalue of the matrix Ra (k), and c^Jmin is the smallest eigenvalue of Rfj (k) 



III. KINETIC GAUSSIAN MODEL ON THREE DIFFERENT FRACTAL GEOMETRIES 



A. The Koch curve, the modified Gaussian model and the master equation 

The fractals that we are going to study are constructed by an iterative procedure in which each segment of 
the object is replaced by a generator. Fig.l shows two different configurations of the Koch curves ||l^ including 
nonbranching, and branching Koch curve. In the iteration, each stage of the iteration is described by a length- 
rescaling factor L, and the number of the segments in the lattice, TV', increases to A'^ by a relation N/N' = L^f , which 
defines the fractal dimensionality Df. Obviously, these examples in Fig.l have different Df, but their topological 
dimensionality Dt is equal to 1 which means that they are quasilinear fractals.. Added to this, another parameter 
which is used to characterize the topological properties of the fractal is R, the order of ramification. The maximum 
and minimum values of i? of a fractal obey the inequality, i?max > 2i?„iin — 2 . 

The examples of the Dt = 1 fractals shown in Fig.l have finite R. The nonbranching Koch curve (NBKC) which 
has Df — In 4/ In 3 and i?min — ^max — 2 shown in Fig. 1(a) is a homogeneous and wiggling chain, while the branching 
Koch curve (BKC) which has Df — In 5/ In 3 and i?min — 2 and i?max — 3 shown in Fig. 1(b) is a inhomogeneous one. 

We assume that the Gaussian spin system with a reduced Hamiltonian 

= fc^cr.o-j, (3.1.1) 

located on these fractals, where /3 = \/kpT,k = J/kpT and the summation is taken over nearest neighbors. 

Unlike Ising spin system, the spin of the Gaussian model can take any real value between (— oo,+oo), and the 
Gaussian-type distribution finding a given spin between Ui and Ui + dui 

f {a,) da, ^ exp [-^'^f^ ^'^^ (3.1.2) 

is assumed to prevent all spins from tending to infinity, where qt is the coordination number of the site i, and bq. is a 
distribution constant independent of temperature. Considering the inhomogeneity of the branching Koch curve, we 
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have assumed that the Gaussian type distribution constants depend on coordination numbers and satisfy a certain 
relation 



(3.1.3) 



As far as we know, this modified Gaussian model appeared firstly in Ref. ||2^ which studied the static critical behavior 
of inhomogeneous fractal lattices. 

In this case the spin transition probability can be expressed as 



where the normalized factor Qi can be determined as 



exp 



i-\-w 



= ^ exp 



exp 



f (ctj) d&i = exp 



and the another useful combination formula can also be obtained 



/oo 
_ -OO 



)da^ 



It." 



(3.1.4) 



(3.1.5) 



So, for magnetic-like perturbation, the master equation suitable for modified Gaussian model on homogeneous and 
inhomogeneous fractal lattices can be written as 



i i \ w / 



(3.1.6) 



B. Nonbranching Koch curve 



First let us focus on the homogeneous nonbranching Koch curve (NB KC) i n which the Gaussian spins are placed 
on the all of sites. Because hq. (t) — h{t) , hq. = b, the master equation ( 3.1.6| ) takes the following form 



= -h it) (^l^^^Yl {^'^1 + ^2 + ^3 + < + Peq {k, {a}) , 



(3.2.1) 



where the mark a denotes generator of NBKC which is shown in Fig. 2(a), the sum goes over all generators, and 
Peg {k, {f}) is the equilibrium distribution function which can be written as 



Peg {k, {<t}) = — exp 



k J2 ''''^1 ^oJ2'^^ 



= 7llexp|/c(aiCT2 +<J2<^3 +'^3<^4: +<^A'^5) 



(3.2.2) 



In Eqs. ( |3.2.lD and ( ^.2.2| ), the coeflicient 1/2 comes from the fact that two neighboring generators share the same 
site 1 and 5. 
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The space-rescaling procedure (see Fig. 3(a)) is performed through the decimation renormahzation transformation 

T" a) = <5 (a*? - O S - a^) . (3.2.3) 

In which the spins ti" , and the interaction k in the a-th generator are replaced by rescahng spin /i'j" , ^2 ^-nd new 
interaction k' , respectively, while the other spins erf 1 '^d, ^'^'^ "'4 ^'"^ integrated from —00 to +00. Under this process, 
the form of the distribution function P^q is invariant. The details of RG calculation is emplaced in appendix (^). 
Here, we give a renormalized master equation 



d\ 1 b{b + 2k) 



dtj e(fc) 62-2fc2 ' ^ \2'^' 2 
1 6^ - Ak^ 



C (k) P - 2fc2 

where 



, ,,,, /&4-4fc2fe2+2fc4\l/2 



k% 

54 _ 4^252 + 2fc4 



(3.2.6) 



Obviously, if the summation for a is arranged in next stage of iteration, the Eq. (3.2.4) can be rewritten as 



1 b{b + 2k) 
dt ) f(kj 62 _ 2fc2 ' 

1 62_4fc2 



h it) E f + + + + J^f ) (fc', {m'}) 



C (fc) 62 - 2fc2 



MO E ( ^/^f + ^^2" + + Mf + ^A^f ) ik', {/i'}) 



(3.2.7) 



Futhcrmore, if we let 



then by time rescaling 



1 6 (6 + 2k) 
" Sjk) 52 _ 2fc2 ' 



1 62 - 4fc2 



C (fc) 62 - 2fc2 1 - 2fel ' 



t _ b^ - ik^b^ + 2fc4 
A/w ~ (6 + 2fc)63 



t = L-H, (L = 3) 



and dynamic parameter transformation 



h (t) ^ h' {t') ^ Xh (t) 



the invariant form of the master equation ( 3.2.2 ) can be restored 

(^')) E {r' + ^2" + + Mf + ) (fc', {m'd 

= (f) f 1 - ^) E (U' + + + A^f + ^A^f ) Pe', ik', W}) 



(3.2.8) 



(3.2.9) 



(3.2.10) 



(3.2.11) 



Let system is in its critical point kc = 6/2 which det ermin ed by the recursion relationship ( p. 2. 6 ), then we can 
obtain the dynamical critical exponent z by use of Eqs. (3.2.9) as 



1 



In 



(6 + 2fc)63 



InL 54 _ 4^^252 _,_ 2fc4 



k^=b/2,L=3 



n4 

= 2— = 2Df. 
In3 ^ 



(3.2.12) 
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However, because 



1 _ 1 Uk' 
V In L V dk 



ka=b/2,L=3 



1 ^-Ak%^-h^^ +2ey 

^ ^ \ - AkW + 2fc4)2 



ka=b/2,L=Z 



In 16 Jn4 



(3.2.13) 



then 



2Df = 2.5237. 



(3.2.14) 



C. Branching Koch curve 



Now, we turn to focus on the branching Koch curve (BKC) which is one of the inhomogeneous fractals. In this case, 
due to the coordination number depends on the place of site, we must to assume that the Gaussian type distribution 
constants depend on coordination numbers and satisfy a certain relation (3.1.3), otherwise, the problem can not be 
solved by applying the decimation RG method directly pl|| . 

In the following we deal with the magnetic-like perturbation master equation that suits for modified Gaussian 
model on inhomogeneous fractal lattices ( [3.1.(: ). We should notice that, for the branching Koch curve (BKC) there 
are two kinds of typical generators (such as a-th and /3-th in Fig. 2(b)): (1) qi = q2 — m — 3, qz = — 2; (2) 
9i = 92 = 94 = 95 = 3, 93 = 2 . For case (1) or case (2), the decimation renormalization group procedure is shown 
in figure 3(b), in which some spins such as cr2, (73 and (74 are integrated, the remainders are rescaled as ^'i and /i2, 
and, at the same time, the interaction k is replaced by k' . Under this process, the form of the distribution function is 
invariant, and the RG transformation of a-th generator is equivalent to /3-th. It can be realized via the calculation of 
the appendix (B 1). 

Our purposes is the renormalization of the master equation. In fact, we only need discuss a typi cal ge nerator. 
Without loss of generality, we take the a-th generator, for instance. The lift and right sides of Eq.( 3.1.6| ) can be 
written as, respectively 



h(t) •A((7)Pe, {k,{a}) 



hs (t) 



h2 (t) 



Peg (fc,{a}) 



(3.3.1) 



hit)-n{k,a)Peg {k,{a}) 
1 . k 

03 



+h3 



3-1 - --2 



+ h3 


^2 - 


- 7-{<^l + f^3 
O3 


^0 


+ /l2 


+ ^3^4 


^5) 


+ h2 


2 ^ 


02 





'^3-7- \P2 + ) 
O2 



(3.3.2) 



where, the coefficient 1/3 (or 1/2) in the terms (or 0-5 ) comes from the fact that three (or two) neighboring 
generators share the same site 1 (or 5 ). 

Multiplying Eqs. ( 3.3.1 ) and ( 3.3.2 ) by the transformation operator 



r(M,f7)-[]<5(/i?-0<5(/i2 -^5), 



(3.3.3) 



and integrate over the {u}, one can obtain (see appendix (B2)) 

i?{h(t) •A(a)Pe, (fc,W)} 
1 



\ la 
Ml 



h'A^) /^^W) TX i^;({fc',M'}) 



M2 



h' (t,fc')-A' {^,k')P' (fc',{/i'}), 



(3.3.4) 
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R{h{t)-n{k,a)Peq (fc,M)} 
= h' {t, k') ■ [R^i (k) . Rn (k)] n' {p', k') (fc', {^'D , (3.3.5) 

in which 

( h'^ {t) h', {t))^{ it) h2 (t) ) (i?A) , 



C b3b2+2kb2-2k^ 2fcb2 \ 

b362-fet'2-2fe2 h3b2-fcb2-2/c2 1 
3fc^ &2(63-fc) I 

/ 4 2 \ 

= ( 3 2 ) , eigenvalues : 5, 1, (3.3.6) 



. _ / (9fc3-l6fc2 62-8/c3)62 ^ 

(.^Ojfc^fc^ - I ^ 9fct-28fe^bi-8fc^62+8fe'' 



(9fc^-16fc2f,2-8fc3)b2 



3 

I 



\ . , 3 3 



, eigenvalues :-, (3.3.7) 



where kc = &2/2 = ^3/3 is determined by the fixed-point equation k* = k' = k { see appendix (B 1)). 
In this case we have to look for the invariant form of the master equation 

{t) A (a) Peg (fc, W) = -h (t) n (k, a) P,^ (fc, {<j}) (3.3.8) 

at the Umit of the order n — s- cx) of the RG transformation. Because our starting point is very close to the fixed point kc 
of the static RG transformation, the eigenvalues of the transformation matrices Ra [k — > kc) and Rfj (fc — > kc) control 
the scaling properties of the largest relaxation time. Hence, according to the foregoing discussion, the dynamical critical 
exponent z is obtained 

_ln(A_/a;„i„)_ln40/3^ 2.3578. (3.3.9) 



Because of 



Ini In 3 



dk'\ _d ( 462^^^ + 62) \ _ 40 



J k=k, dk V 8fc4 - 8fc362 - 286p2 + 954 J ^^^^ 3 



1 _ 1 Uk' 
V InL V dk 



k=ka 



In 40/3 
ln3 ' 



(3.3.10) 



then 



z = - = 2.3578. 

V 



(3.3.11) 



D. Multibranching Koch Curve 

Based on the branching Koch curve, we now construct another generalized one, the multibranching Koch curve 
(MBKC), and investigate its critical dynamical behavior of the kinetic modified Gaussian model on this lattice. The 
constructional process is shown in Fig. 4. Obviously, it also is a inhomogeneous example of the Dt — 1 fractals which 
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has Df = ln(2m + 3)/ln3, i?min = 2 and i?max = m + 2, m = 1, ...oo. The effect of the order of ramification, R, on 

the critical slowing down could be seen by this example. 

Similarly, there are two kinds of typical generators in the multibranching Koch curve: (1) qi = q-s = q^ ~ l/(m + 2), 
= 2, (j = 5, • • • , m + 4), 92 = 2 ( a-th generator); (2) qi = 53 = 54 = (j2 = l/(m + 2), = 2, (i = 5, • • • , m + 4) ( 

f3—th generator). In fact, the decimation renormalizing procedures of t hese two cases result in the same consequence 

(see Fig. 5). It can be realized via the calculation of the appendix (CI). 
Our purposes is the renormalization of the master equation 



d_ 

dt 



(3.4.1) 



or 



dt 



h (t) -A (a) Peg (fc, {a}) = -h (t) n (fc, a) P^g (fc, {a}) . 



(3.4.2) 



In fact, we only need discu ss a typical generator. Without loss of generality, we take case (1) for instance. The lift 
and right sides of Eq.( 3.4.1 ) can be written as, respectively 



hit) -A (a) Peg {k,{a}) 



l'm+2 



1 



TO + 2 



„OL I „a I r 



^ m+4 



. i=5 



Peg (k, W) 



(3.4.3) 



and 



h{t)-n{k,a)Peg (fc,M) 



'yl l'm+2 

+ hrn+2 
m+4 



TO + 2 ^ 6, ^ 



Jm+2 



l'm+2 



m+4 



bm+2 



m+4 



4 ~ 2^^^ + '^2 



Wm+2 



i=5 



02 



+ /l2 



^eg (fc, W) . 



(3.4.4) 



By virtue of ( C.2.1 - C.2.5 ), the results of the decimation RG transformation of ( 3.4.3 ) and ( 3.4.4 ) are, respectively 



i?{h(i) ■A{a)Peg (fc,W)} 
1 

7! 



{h'^^,{t,k') K[t,k'))[ r+f" )p;(fc',{A*'}) 



2M2 



\^{t,k')-M {^,\k')p'{k',{^,'}) 



(3.4.5) 



where 



( h'^+2 K) = { hm+2 h2 ) (i?A) 



/ b2b,Ti + 2-62fc-2mfc-'+fcb2(™+2) 262 fc 

(-/^a) I {b2bm+2—b2k—2mk^) b2bm + 2—b2k—2mk^ 



{m+2)mk'' 



62 



-62^— 2mfc^ 



fc2 6„ 



b2bm^2—b2k 

-b2k—2mk^ 



and 



i?{h(i)-f2(fc,f7)Peg (fc,W)} 

h' «, f) ■ [RJ' (k) ■ Rn (k)] ■!!' (,.', f) P.', (k', {,.'}) 



a/ ,,fa 

P2 



fcm + 2 

2/^2 



P'egik'AP-'}) 



(3.4.6) 
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where 



^ 2k' 
{n 



1 / I ' V y V 021 a22 



_ -2bm+2mk^ - 6„i+2&2fc^ - + ^2&to+2 - &2fc^ (to + 2) &„i+2 + fc'' (m + 2) m 

" ib2bm+2 - hk - 2mk'^) bjn+2 (k + b„i+2) 

2fc3 (mfc + 62) 



ai2 



021 



&m+2 {k + b.m+2) {2mk^ + b2k - 62V+2) ' 
{m + 2) {mk + 62) 



022 



62 (A; + bm+2) {2mk'^ + b2k - 52&m+2) ' 

262mfc-'' + b^P + 262&m+2mfc^ - bi6^„+2 + 2b„,+2&2fc^ - 2fc^m 
^2 (fc + bm+2) (2mfc2 + feafc - b2bm+2) 



\ . . m + 2 m + 2 



= ( ^'7'^ ^ ) ' ^ ^(^' ^(^^TT^ (3.4.7) 
/ m + 3 2 \ 

{Ra)i.^i. —{if , n\ ,1 , ejoenwaZues : 2m + 3, 1 (3.4.8) 

/fc^fce \^ im (to + 2) TO+ly ^ ' 

where fee is determmed by the fixed-point equation k* — k' — k (see ( |C.1.12 )) as 



&^ ^ W2 ^ ( ^ 

2 TO-t-2 ^ ' 



From Eqs. (3.4.5) and (3.4.6) we can see that, in order to look for the invariant form of the master equation 



jh (t) A (a) Pe, {k, {a}) = h (t) U (fc, a) {k, {a}) (3.4.10) 

we need to do the renormahzation endlessly up to the limit of the order n ^ cxd of the RG transformation. But, 
because the system is very close to the fixed point kc of the static RG transformation, the scaling properties of the 
largest relaxation time are under the control of the largest eigenvalue Amax — 2to + 3 of the matrix Ra {k — > kc) 
and the smallest eigenvalue t^min = 2(m+3) ^ ^'^^ • Hence, the invariant form of the master equation can be 

obtained through preforming the time rescaling 

t^t' ^ L-H= , (3.4.11) 

'^max/ ^min 

and from here, the dynamical critical exponent z can be got 

_ In (Amax/ l^min ) _ 1 2 (2to + 3) (to + 3) ^ 

h^L ^ hTs TO + 2 ■ y-- ) 



Because of 



dk' 
Ik 



462fc3 (&2 + mk) 



= 2 



* \8fc4TO - 862m/c3 - 4 (to + 2) blmk"^ - ABbl {m + 2)- Ak^bl + (to + 2)^ fo| / ^^^^^^ 
(m + 3) (2to + 3) 



TO + 2 
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then 



In 3 



■In 



2(m + 3) (2m + 3) 



(3.4.14) 



IV. CONCLUSIONS 



Based on the dynamical real-space renormalization proposed by Achiam and Kosterlitz, we have suggested a gen- 
eralizing formulation that suits arbitrary spin systems. The new version replaces the single-spin flipping Glauber 
dynamics with the single-spin transition dynamics. As an application, we focused our minds on the kinetic Gaussian 
model (cr = — oo, ...,oo, continuous spin model), and studied three different fractal geometries with quasilinear lattices, 
including nonbranching, branching and multibranching Koch curve. We calculated the dynamical critical exponent 
z for these lattices using an exact decimation renormalization transformation in the assumption of the magnetic- 
like perturbation, and found that it can be written universally as z = ^/v, where v is the static length-correlation 
exponent. 

In the first example, the nonbranching Koch curve, z — 1/ v — 2Df, Df — In 4/ In 3 is the fractal dimensionality 
of the NBKC. Being a quasilinear chain, the geometrical effect of the wiggliness of the nonbranching Koch curve is 
that the correlation length C of the one-dimensional linear chain should be replaced by the real correlation length ^, 
( = (V^/ However, for one-dimensional linear chain with Gaussian spin each lattices, we have known that the 
critical dynamical exponent 2 = 20. Sot~C^ = (C'^''^ )^ = means z = zDf — 2Df. This result is coincided with 
our calculation by DRSRG technique. 

In the branching Koch curve, the result z — l/i/ is also valid, but that i = In ^ = ^'^^^^ is not simply 

related to the fractal dimensionality Df. 

In the multibranching Koch curve, the result z = 1/j/ is obtained once again. We can see from z — 
j-^ In ^("'+'^)(^™+'^) that, when m ~ 0, z = 2, the lattice is one-dimensional chain and z is equal to the rigorous 

result; when m = 1, z = corresponding to the branching Koch curve; when m — > oo, z — > oo (see Fig. 6). All 

of these means that the critical slowing down of the Gaussian spin on Koch curve is dependent mightily on the order 
of ramification R. 

In fact, in our previous paper we have found that for a translational symmetric lattices with Gaussian spin model, 
the critical dynamical exponent z = 1/v, v = 1/2 at the critical point Kc = b/2d based on rigorous calculation [Q. 
Yet, in this paper the result z = l/i/, have been proved once again by the dialational symmetric lattice systems. 
We guess that z = could be a universal conclusion for kinetic Gaussian model. Of course, we must realize that 
the result what we have obtained in this paper is carried out in the assumption of the magnetic-like perturbation. 
However, the perturbation itself (magnetic-like or energy-like) is only a special assumption. For a general perturbation 
the master equation is not always invariant under the RG transformation because the perturbations probably have 
components along all the relevant operators. By this token, whether the z — 1/v will be a universal conclusion for 
kinetic Gaussian model waits for further investigation. 
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APPENDIX A: RG CALCULATION OF NBKC 



To perform the decimation transformation, we have only to multiply both sides of the master equation ( |3.2.1 ) by 
the transformation operator 



11 



and integrate over the {a}, i.e. 

(1) M^) E ^ { i^r'' + + + + ^^^^ ^""^^l 



where 



RiafPe, (fc,W)} 

r daidc72 ...dcTN^lsUl- 4) 6 U - 4) afPeg {k, {a}) 

— 00 o 



= = R{Pe, {k,{a})}-^ 



(A.1) 



(A.2) 



(A.3) 



/OO 
da^da^da^ exp {A; (/U^ct^ + a^a^ + a^a^ + </U2 ) 
-00 

' ' f k^ 



6(62-2fc2)^''P 1 6(62-2fc2) 



1 (5- 



2 /l2 _ i,2\ 



2 6 (62 - 2fc2) 



/OO 
da^da^da^ ■ exp {fc (^^j'cr^ + ct^ct^" + alj'crj + trj /i^) 
-00 



(^2f + K)' + Kr 



62 - 2A;2 



fc2 



62 - 2fc2 



M2 



fc2 



62 - 2fc2 



W. 



The remanent integration i?{Peg ik, {a})} is a important one. We hope to keep on a invariant form of the trans- 



formational distribution function P'. 



^ ^ / '^'^l'^'^2 • • • rfcTjv n ^ ('"1 ~ '^l ) ('"2 - ) exp 
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P 
4 



1 / 



o Ml 



(27r) 



■ exp 



^ y V ^ (62 - 2fc2) ^ [ 6 (62 - 2fc2) 

1 / «\2 1 



6 6* - 4fc262 + 2k^ 

'2 62(52_2fc2) 



("0 



/9 /3 



Obviously, one must to rescale the spins /i", /if and interaction k so as to keep the equihbrium distribution function 
to be invariant 



52(62_2fc2) ) 



(A.4) 



fc' 



fc36 



54 _ 4^252 + 2fc4 



then 



(A.5) 



i?{Pe,(fc,w)} 



(27r)-^ 



) 1 I IP 10 

ZllV6(62-2fc2)^"Pi'^^i ^2 "2 

= exp i fc' ^ \ (y-i) 



2 1 / 2- 



2lM2 



.2 1 / ,^^2- 
Ai2 



(A.6) 



Eq. (A.5) is reputed to be the recursion relation which enable one to determine the fixed point of the static RG 
transformation k^. 
Upon that, we have 



i?{arPe,(W)}-/i?Fl (fc',M) 



(A.7) 



R{aqP,, (fc,M)} 



Mr 



62 - 2fc2 



(A.8) 
(A.9) 



R {a-P,, (k, {a})} = (/.? + M^) P,', (fc', {/i'}) , 



PK"Pe, (fc,W)} = 



fc2 



/^2 



62 - 2A;2 



K+M2) 



Substituting Eqs. ( [a1| - |A.11| ) into ( |A^ ), one can obtain 



d\ 1 b{b + 2k) 



dt J ^{k) -2k^ ^ \2''' '2 

1 62-4fc2 
'C(fc) 62-2fc2'^'"'^ V2 



(A.IO) 
(A.ll) 



(A.12) 
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Obviously, if the summation for a is arranged in next stage of iteration, the Eq. (A. 12) can be written as 
d\ 1 b{b + 2k) 



e (fc) 62 - 2/c2 
It is just the Eq. ( ^.2.'^ ). 



(A.13) 



APPENDIX B: RG CALCULATION OF BKC 



1. The RG transformation of a-th generator is equivalent to /3-th 

We can show that the RG transformation of a-th generator is equivalent to /3-th, but the precondition is that the 
Gaussian type distribution constants depend on the coordination number and satisfy the relation (3.1.3). It can be 
realized via the following calculations. 

The effective Hamiltonian of the /3-th generator is 



1 



KbT 



{a, k) — k {^i(J2 + + 



13 , 13 , 13 



4"5 



63 

2 



3 



(B.1.1) 



where, the coefficient 1/3 in the terms {^i^ ^-nd ^"'5^ comes from the fact that three neighboring generators share 
the same site 1 and 5. We take the decimation renormalization transformation operator as 

then, by integrating spins (T2, fa and (T4 from — cx) to -fcx), one obtain 

1 



where 



R lexp 



/oo 
d(Jida2 . . . da^ • a) exp 

-00 



C exp \ k^^il^ - y ^ 



c = 



h^hl (1 - - 2fc2/6263) ' 



(B.1.2) 



(B.1.3) 



k^{k + b2) 



(6263 - kh2 - 2fc2) (k Jr 63) 



if taking 



3fc* - 2k^h3 - 4fc26263 - 26^fc2 + 6263 

hihh^kbr^WUk + bsj ' 



/- _ /3fc4-2Pfe3-4fc26263-262fc2+6253 
^ " V 63 (M3-A:62-2A;2)(fc + 63) ^' 



(B.1.4) 
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ko _ fc3 (fc + bz) ^3 

C ^ 3fc4 - 263/23 - Ak^b2b3 - 2k%l + b2bl ' 



(B.1.5) 



then 



R |exp 



1 



(B.1.6) 



The same as above, the effective Hamiltonian of the a-th generator is 

^e/ / K) = fc (O-^ (72 + ^3 + ^2 ^4 + ^^3 ^4 + ^4 ^^5 j 



63 
' 2 



1 



62 
2 



(B.1.7) 



where, the coefficient 1/3 (or 1/2) in the terms (u")^ (or (0-5 )^ ) comes from the fact that three (or two) neighboring 
generators share the same site 1 (or 5 ). By integrating spins (72, CT3 and (74 from —00 to +00, one obtain 

1 



R < exp 



KbT 



Kff{<J,k) 



1 



/OO 

d(7id(72 • . . d<7^ • (fi, a) exp 
'OO 

Cexp jA;oM?/i? - (a^?)' - {f^^f 



(B.l. 



where 



fc3 (fc + 62) 



(6263 - kb2 - 2fc2) (fc + 63) ^ 



a = 1 - 3fcV&3 



6 = 1- 



3fcV^ 3fcVb26i 

1 - fc2/6i " (1 - fc/53) (1 - fc/53 - 2fc2/fe263) ' 

2kVblbl 



2kyb2b3 

1 - kybl (1 - fc/63 - 2fc2/6263) (1 - k/bs) ■ 



Based on the relation (3.1.2), i.e.. 



we can see 



So, if only 



Ci = 6 = e = 



62/63 = 2/3, 



8fc4 - 8fc362 - 286p2 + 954 
62 (352 - 2fc62 - 4fc2) (2fc + 362) 



/ 8fc4 - 8fc362 - 2862fc2 9^4 
62 (362 _ 2kb2 - 4fc2) (2fc + 362) ^' 



(B.1.9) 



fc' 



fco _ 462fc^ (fc + 62) 

e ~ 8fc4 - 8fc362 - 286ifc2 + 9bi ' 



(B.l. 10) 



one can get 



R 



ia,k) 



CexpjfcV'i^^f 



62 1 
2 2 



(B.l.ll) 
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In fact, ( |B.1.4| - pX^ ) is coincided with ( B.1.9 - B.1.10 ). Solving fixed-point equation k* ~ k' = k, the 
is obtained 

^ _ &2 _ ^ 

Now, as the decimation renormalization transformation operator is taken as 



from (B.1.6) and (B.1.11) we can see 



i?{Pe, (fc,M)} = -exp 



PL {k',{^i'}) 



This means that the distribution function is invariant under RG transformation. 



2. RG transformation of master equation 



Besides R{Peq (fc,{(T})}, one can also calculate 



R {a^P,, (fc, {a})} = E {T (/i, a) a^Pe, {k, {a})} 

-Pe, (fc',{M'})--^ 



where 



/oo 
da^da^da^ exp {fc (^^fi^ + a^a^ + cr^^a^ + a^af + a'^ fi^) 
-OO 



-OO 

' 2 



exp 



■ (6263 - ^62 - 2k^) (63 + fc) 



fc3 (62 + A:) 



(6263 - kb2 - 2fc2) (63 + k) 



M>2 



da^d(7<^d(72 ■ af exp {fc (^i^CT^ + ct^ct^ + cr>^ + ct^< + crj/i^) 



2 



y('X3) 



Upon that, we have 



R {a?Pe, ({fc, a})} = ^^r^; {{k', : 



R{0^Pe, ({fc, O})} = -^f^TPL {{k', m'}) : 



"P ^/i. 1 fe(M2-fc2X + fc^(fc + fe2)/i^" „, 
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RW^Pe, i{k,a})} 



^b2b3-kb2-2k^ 



By virtue of these integral results, Eqs. ( |3.3.4 ) and ( 3.3.5 ) can be obtained. 



(B.2.7) 
(B.2.8) 



APPENDIX C: RG CALCULATION OF MBKC 



1. The RG transformation of a-th generator is equivalent to /3-th 

For MBKC, we can also show that the RG transformation of a-th generator is equivalent to /3-th, but the precon- 
dition is, the same as about BKC, that the Gaussian type distribution constants depend on the coordination number 
and satisfy the relation (3.1.3). It can be realized via the following calculations. 



The effective Hamiltonian of the a-th generator is 
1 



m+4 



CTl CTg + CTg (74 + (T4 2^ ((74 + (73 ) (7, 



j=5 



m -I- 2 



&2 

2 



'm+4 



i=5 



(C.1.1) 



where, the coefficient 1/ (m + 2) (or 1/2) in the terms (crf)^ (or (ctj)^ ) comes from the fact that (m + 2) (or 2 ) 
neighboring generators share the same site 1 (or 2 ) . We take the decimation renormalization transformation operator 
as 



T" (A.,(7")=,5(/.?-(7n<5(A^^-(72"), 

then, by integrating spins (73, (74 and Oi (« = 5, • • • , m + 4) from —00 to +00, one obtain 

1 



where 



R lexp 



KbT 



/oo "1+4 
n d(7f -T" (M",a")exp 

C exp < 



1 



KbT 



Kff 



a,, a i>m+2 1 ^ / „x2 ^2 1 ^ . „^2 



)ri/2 



(27r)' 62 



{brn+2 + fc) (&2&m+2 - fefc - 2mfc2) 



(C.1.2) 



(C.1.3) 



fco = 



k'^ (62 + mk) 



{bm+2 + fc) (62fom+2 - ^2^ - 2mfc2) ' 



m{m + 2)k^ - 2bm+2mP - 26^+2"^^^ - ^2^:^ (m -|- 2) 6,„+2 - fc^fe2&m+2 + &m+2^2 

^m+2 (6m+2 + (fc2&m+2 - ^2^ - 2TOfc2) 



6 = 



-2fc2&2&m+2 + 2fc''m - 2b2bm+2ink'^ + &m+2^2 - 2b2mk^ - 

^2 (&m-|-2 + ^) (&2&m+2 - ^2^ ~ 2mfc2) 
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Base on the relation (3.1.3), i.e 



we can see 



bm+2 m + 2 



e = Ci = 6 

-Sfc^^TO + 8b2mk^ + 4 (m + 2) blmP + APbj (m + 2) + APbj - (m + 2)^ b 



2f,4 
2 



62 ((m + 2) 62 + 2k) {-bl (m + 2) + 262^ + Amk^) 



(C.1.4) 



(C.1.5) 



4fc^ (&2 + mfc) 



if taking 



then 



((m + 2) 62 + 2fc) (6^ (m + 2) - 2^2^ - Amk"^) ' 



' fc I,' 



i? I^exp 



2 m + 2 ' 2 2 
The same as above for case (2), the effective Hamiltonian of the (3-th. generator is 
1 



m+4 



" l,(.f)%(.f)V(4)V;;^(.0' 



Wm+2 



62 

2 



■m+4 



(c.1.6) 



(C.1.7) 



(C.1.8) 



(C.1.9) 



where the coefficients 1/ (m + 2) in the terms (^^i^ ^-nd (^"'2) come from the fact that m + 2 neighboring generators 

share th e same site 1 and 5. By integrating spins 0-3, a a and Gi (« = 5, • • • , to + 4) from —00 to +00 and using the 
relation (C.1.4), the same result can be obtain 



R ^exp 



KbT 



/oo "i+4 
nd'xf .r^(M,a^)exp 
-c>o „-_i 



C exp < k Hi IJ.2 



I ip ip "m+2 



1 



TO + 2 



1 



TO + 2 



'/3 
M2 



Therefore, from (C.1.8) and (C.l.lC), we have 



i?{Peg (/c,W)}--exp 



where the recursion relation is 
fc' = - 



4fc2 (^2 + TOfc) 62 



4to - 862m/fc3 - 4 (to2 + 3m + 3) bjk^ + (to + 2)^ b^ 



(C.1.10) 



(C.1.11) 



(C.1.12) 



Expression (C.1.11) means that the distribution function is invariant under RG transformation. 
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2. RG transformation of master equation 



Besides R{Peq (k, {c})}, one can also calculate 

R {afP,, {k, {a})} = ^ {T (m, a) a^P^, {k, {a})} 



/oo 
daida2 . . . d(7Ar n '5 (m? - O ^ (/x^ - a^^) afPeg {k, {a}) 
-OO a 



w 



where 



W 



POO /m+4 \ f 

J da^da^ f J| dcrfj exp Ik 



m+4 



K)' + K)' 



62 
2 



27r 



1 

+ - 



m/2 



k'^ (62 + '77,A:) 



■?n+4 



i=5 



^ (fem+2 + fc) (&2&m+2 " ^'2^ " 2mF) \ (6^+2 + A:) (62^m+2 " b^k - 2mk'^) 
k^ (62 W2 - mk"^) 



2 (&m+2 + A;) (62&m+2 - ^2^ - 2mfc2) 



bm+2 



J— 5,...,m+4 



62 
2 



i=5 
■m+4 



m+4 



Ml 0^3 + 0-3 0^4 + Cr4 M2 + 2^ (0^4 + 0^3 J 



j=5 



L i=5 



/oo / m+4 ^ „ 



dcr" ■ cr" exp 



da" exp 



M< + ^3)^,"-|«)' 



X exp <; fc + c^3 + ctJ/x^] 

fc2 



Om+2 



b2bm+2 — b2k — 2mk'^ 



W, 



/oo /m+4 
-°° V i=5 



CT4 exp < k 



Ml 



m+4 



CT3 + CT3 CT4 + /Z2 + 2^ (0-4 + CT3 ) (7, 



l^m+2 



&2 

2 



'm+4 



E 



j=5 



(m/c^ + 62^2) + (fc2&m+2 - ?T^fc2) fc/i^ 

(/c + bm+2) {b2bm+2 - b2k - 2mk'^) 



i=5 



/oo /m+4 \ ( 

daedal ( []'^<l'^3exp|fc 



m4-4 



Ml 0^3 + 0^3 0^4 + 0-4 M2 + 2^ (<^4 + 0^3 ) 



i=5 
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and then 



Wm+2 



2 



'm+4 

. i=5 



{{mP + b2k) + {b2b,n+2 - mk^) fi'^) k 
{k + bm+2) (&2fcm+2 - hk - 2mk^) 



R {a?p,, (fc, M)} = ^^^'^Pe, (fc', {^^'}) 



R{a^Pe, {k,{a})} = ■^^^TPeg ik',W}) 



(C.2.1) 
(C.2.2) 



{k + bjn+2) (020m+2 - 02k - 2mk^) 



(C.2.3) 



V'f + Om+2) [b2bm+2 - 02fc - 2mfc^) 



(C.2.4) 



i?{a^"Pe,(fc,W)} 



j=5,...,m+i 1 



V? ^2bm+2 -b2k - 2mk'' 



(C.2.5) 



By virtue of these integral results, Eqs. (3.4.5) and (3.4.6) can be obtained. 
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FIG. 1. Three stages in the construction of fractals. The generator of the fractal appears in the first line. The second and 
the third lines correspond to the second and thirds stages of iteration. Different fractals are placed in different columns: (a) A 
Non-branching Koch curve (NBKC) with Df = In 4/ In 3 and R = 2. (b) A branching Koch curve (BKC) with D/ = In 5/ In 3, 

Rmin = 2 and -/?max = 3. 
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FIG. 2. The second construction stage of the NBKC and the BKC. (a) NBKC: All of the generations are the same; (b)BKC: 
there are two kinds of typical generators (such as a-th and /9-th). 

FIG. 3. Decimation RG procedure :(a) NBKC; (b) BKC. 

FIG. 4. The construction procedure of multi-branching Koch curve (MBKC) with Df = ln(2m -|- 3)/ In 3. 

FIG. 5. Decimation RG procedure of the multi- branching Koch curve (MBKC). 

FIG. 6. Critical dynamical exponent of the kinetic Gaussian model on the Multi-branching Koch curve. 
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